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We introduce the concept of distillability sudden death, i.e., free entangled states can evolve into 
non-distillable (bound entangled or separable) states in finite time under local noise. We describe 
the phenomenon through a specific model of local dephasing noise and compare the behavior of 
states in terms of the Bures fidelity. Then we propose a few methods to avoid distillability sudden 
death of states under (general) local dephasing noise, so that free entangled states can be robust 
against decoherence. Moreover, we find that bound entangled states are unstable in the limit of 
infinite time. 
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I. INTRODUCTION 

Entanglement is not only a remarkable feature which 
distinguishes the quantum world from the classical one 
but also a key resource to realize high-speed quantum 
computation and high-security quantum communication 
In realistic quantum-information processing, entan- 
glement usually needs to be prepared or distributed be- 
forehand among different remote locations. However, in 
the process of entanglement distribution the quantum 
systems are not isolated and each system will unavoid- 
ably interact with the environment. This leads to local 
decoherence which will degrade the entanglement of the 
shared states. It is thus of fundamental importance to 
study the entanglement properties under the influence of 
the local decoherence. In this context, Yu and Eberly 
[2| investigated the time evolution of entanglement of a 
bipartite qubit system undergoing various modes of de- 
coherence. Remarkably, they found that, although it 
takes infinite time to complete the decoherence locally, 
the global entanglement may vanish in finite time. The 
phenomenon of finite-time disentanglement, also named 
entanglement sudden death(ESD), unveils a fundamental 
difference between the global behavior of an entangled 
system and the local behavior of its constituents under 
the effect of local decoherence. Clearly, ESD puts an limit 
on the applicability of entangled states in the practical 
quantum information processing. 

Initially, Yu and Eberly reported the ESD for two- 
qubit entangled states, but this effect is not limited to 
such case. Further investigations in a wider context 
including higher dimensional Hilbert spaces have been 
made byvarious groups 0, i, H, i, 0, S, i, M, [HI, E, M, 
[3,[ll,E3, 17, 18, i^. There are also a number of studies 
looking at ESD in more complicated systems using other 
entanglement measures [13, HU, [H, [H, [H, III, [l^l. 
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and an attempt to give a geometric interpretation of the 
phenomenon has also been made [2^. In addition, exper- 
imental evidences of ESD have been reported for optical 
setups and atomic ensembles (30l |. 

However, all previous studies omitted an important 
fact that high dimensional bipartite entangled states can 
be divided into two classes [3l|. One is free, which 
means that the state can be distilled under local oper- 
ations and classical communication (LOCC); the other 
is bound, which means that no LOCC strategy is able 
to extract pure-state entanglement from the state even 
if many copies are available. Bound entanglement (BE) 
cannot be used alone for quantum information process- 
ing, and irreversibility occurs in asymptotic manipula- 
tions of entanglement for all BE states [32|- Since it was 
constructed from a pure mathematical point of view, we 
may ask whether BE can appear in physically relevant 
quantum systems naturally. Very recently a few works 
have addressed this question [s^. Their results suggest 
that different many-body models present thermal bound 
entangled states. 

In this paper, we investigate the problem from a very 
different viewpoint, i.e., in the present of local decoher- 
ence, which is one of the dominant noises during the 
distribution of entanglement. Analogous to the defini- 
tion of ESD, if an initial free entangled state becomes 
non-distillable in finite time under the infiuence of local 
decoherence, then we say that it undergoes distillabil- 
ity sudden death(DSD). Note that when a free entangled 
state loses its distillability at a specific time, it may still 
be entangled since we cannot exclude the existence of 
BE. So far, it is not clear whether bipartite bound en- 
tangled states can be created from free entangled states 
under local decoherence process naturally. The first aim 
of this paper is to show that such a process indeed exists 
through an explicit qutrit-qutrit example. 

Afterwards we propose the DSD-free state, which has 
entanglement robust against local decoherence. Such 
entangled states are thus useful resources for practical 
quantum-information processing. So the second aim is 
to address the DSD-free state. We develop a few sys- 
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tematic approaches to build DSD-free states. Finally we 
will show that different from free entangled and separa- 
ble states, no PPT bound entangled states exist in the 
infinite time limit. 

The paper is organized as follows. In Sec. II, the idea 
of DSD is presented. In particular, we show the phe- 
nomenon of DSD through a specific qutrit-qutrit exam- 
ple. In addition, we show that one can avoid the sudden 
death of distillability by performing a simple local unitary 
operation on the initial state. Furthermore, we compare 
the behaviors of states affected by decoherence in terms of 
the Bures fidelity. In Sec. Ill, we develop some methods 
to build DSD-free states which can protect free entangle- 
ment under general local dephasing noise. We also show 
that no PPT entangled states can exist in the limit of 
infinite time. Finally, in Sec. IV we discuss some open 
questions and also give a summary of our results. 

II. QUTRIT-QUTRIT DSD STATES UNDER 
LOCAL DECOHERENCE 

Before discussing dynamical process of entanglement, 
we briefiy review how to characterize bound entangled 
states. It was proven in [3l| that a quantum state with 
a positive partial transposition (PPT) is non-distillable 
under LOCC. Therefore PPT entangled states must be 
bound entangled states [s^l- To verify them, one can use 
the so-called realignment criterion (cross-norm criterion) 
[35| . The definition of realignment on the density matrix 
is given by j,; = Pikji- A separable state p always 

satisfies ||p^|| < 1- For a PPT state p, the positive value 
of the quantity — 1 can hence verify that it is a 

bound entangled state. 

The system we study consists of two noninteracting 
qutrits in two independent local environments, each cou- 
pling to one of the qutrits. Here the qutrit is a three- 
dimensional state composed of the computational basis 
|0), |1), and |2). Besides, we take the weak local de- 
phasing noise to model the environment, such noise is 
indeed one of the main decoherence sources in solid-state 
systems. 

The general time-evolved density matrix expressible 
in the operator-sum decomposition is the completely 
positive trace preserving map [sBl p{t) = e(p(0)) = 
^1 (.t) P (O)-f'^M (0- The operators {K^ (<)} represent- 
ing the infiuence of statistical noise satisfy the complete- 
ness condition (t) Kj^ (t) = I which guarantees 
that the evolution is trace-preserving[l|. In our model, 
the operators are of the form Kfj,{t) = Dj{t)Ei{t) such 
that 
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p{t)^£ {p (0)) = ^ d] (t) El [t] p (0) it) D, (t) .(1) 

Here, E.^ (t) and Dj (t) correspond to local dephasing noise 
components acting on the first and second qutrits, re- 
spectively, and both operators satisfy the completeness 
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FIG. 1: (Color online). The eigenvalue Ai as a function of 
time t and dephasing rate V for two different cases: (a) a = 
4.5; (b)Q! = 4.9. Figures (c) and (d) correspond to the contour 
plots of the upper one, respectively. 



conditions. For simplicity, we first take these to be of the 
specific forms 

Ei{t) = diag(l,7A,7A)«)l3 , 

E2{t) = diag(0,WA,t^A) «il3 , 

Di{t) = I3 (g) diag(l,7B,7B) , 

D2{t) = I3 (g) diag(0,WB,WB) , 

where I3 is the 3x3 identity matrix, 
7A (t^^ e-rA*/2^ 7B(t) = e- ^^-'/^ uo^{t) = 

v/1-7aW, and lob {t) - VT^tIW ■ 

For concreteness, we illustrate our ideas by considering 
the following qutrit-qutrit state, 

P(0) = ^(|01) + |10) + |22))((01| + (10| + (22|) 
+ ^(|00)(00| + |12)(12| + |21)(21|) 
+ ^ (111) (111 + 120) (201 + 102) (02|) (2) 

with 4 < a < 5. In fact, it is straightforward to prove 
that the state p (0) is free entangled. By using the lo- 
cal operation (|0) (0| + |1) (1|)«)(|0) (0| + |1) (1|), one can 
converts p (0) into a 2 2 entangled state, which is af- 
terwards distillable by virtue of the BBPSSW-Horodecki 
protocol [131 • So the state p (0) is a free entangled state 
and we take it as the initial state under the local dephas- 
ing noise in Eq. (1). This evolution can be calculated 
analytically, i.e., at time t the two-qutrit density opera- 
tor p{t) reads 
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{|00),|01),|02),|10) Jll>,|12),|20),|21),|22)}. 

Let us analyze the matrix carefully. There are three 
eigenvalues of the partial transposition of the state p (t) 
which could be negative, namely Ai = / (r^) , A2 = 
/(Fb), and A3 = fiTA+TB), where /(A) = 
%^ (l05e^* - ^11025e2^* + 1764e^* (4 - 5ae^* + a^e^*) 

For simplicity, we choose the local asymptotic dephasing 
rates Fa = Tb = F, and thus Ai = A2 < A3 in the follow- 
ing arguments. In order to have a vivid illustration, in 
Fig. 1 we plot Ai as a function of t and F with specific 
a = 4.5 and a = 4.9, respectively. Fig. 2(a) shows the 
value of Ai versus t for different decoherence rates with 
a = 4.5, and we can see that the eigenvalue of the partial 
transposition of the state p{t) will always arrive at a 
positive value in finite time. For example, if we choose 
F = 1, the density matrix p{t) will become a PPT state 
after time t w 0.58 in Fig. 2(a). Analytically, the time at 
which p{t) becomes a PPT state is td = ^ In ^^^^^^ . 

Next, we use realignment to verify the BE in 
this evolution. To this end, we need to com- 
pute the quantity ||p(t)^|| — 1, and it is given by 
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order to make a comparison, we plot pit)' 
a function of time t by choosing a ~ 4.5 and F 
in Fig. 2(b). We can see that if i < 0.84, the value of 

p{t) — 1 is always positive, which indicates that in 

the range 0.58 < t < 0.84, the two-qutrit system is a 
bound entanglement state. Thus we have shown that 
free initial entangled states can evolve into non-distillable 
(bound)entangled states in finite time under the local ex- 
ternal asymtotic dephasing noise. 

Besides we find that with a = 4.5 and F = 1, the 
state p{t) will be separable when t > 1.39. It could 
be proved by extracting three "2 x 2" density operators, 
which are respectively spanned by {|00) , |01) , |10) , |11)}, 
{|01),|02),|21),|22)} and {|10) , |12) , |20) , |22)}, from 
the state p{t). Then one can find out that the condi- 
tion t > 1.39 makes the three states separable. As the 
state p(t) is actually the linear combination of the above 
three states and some product states, p{t) also becomes 
separable when t > 1.39. Nevertheless, there is still a 
small window t £ [0.84, 1.39] in which the state p{t) has 
merely PPT and we don't know whether it's entangled 
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FIG. 2: (a)The eigenvalue Ai versus time t for different de- 
phasing rates. The solid, dash-dotted, and dashed lines cor- 
respond to r = l,r = 0.7, and F = 0.4, respectively, (b) 

The quantity |jp (*)^|| ^ 1 versus time f for F = 1. The other 
parameter a = 4.5. 



or separable. 

To summarize the example, start from the initial state 
p(0) in 4 < a < 5 the qutrit-qutrit state p{t) will ex- 
perience three phases under the local dephasing noise in 
Eq.(l): it evolves from an initial distillable entangled 
state to a PPT entangled state, and then changes to be a 
separable state and finally remains in this phase for even. 
In other word, p{t) first experiences DSD and then ESD 
under local decoherence. 

In contrast with the above example, we briefiy consider 
the initial state p(0) for 3 < a < 4. In this case, it is a 
PPT entangled state (38l |. one can verify that the state 
p{t) will finally become separable in finite time under 
the local dephasing noise in Eq.(l); that is, p{t) only 
experiences ESD. Moreover, we can easily verify the time- 
domain factorization relation as follows, 



til 



(4) 



where ii,t2 > 0, and S denotes the noise in Eq.(l). Be- 
cause the entanglement cannot increase under the local 
operations, we conclude that if a state becomes separable 
at a specific time in Eq.(l), it must remain separable in 
all subsequent time. More generally, if an entangled state 
evolves into a PPT state (either entangled or separable) 
at a specific time in Eq.(l), it always has PPT for all 
subsequent time since the local operation cannot change 
PPT. 

Of course, such disappearance of distillability in finite 
time can seriously affect the application of entangle- 
ment in quantum information tasks. An important 
question arises naturally: given a free initial entangled 
state, does there exist a suitable intervention that may 
alter the final fate of DSD? In the following we show 
that one can realize this aim by merely performing 
a simple local unitary operation on the initial state 
in the above example. Let the initial local operation 
he U = \, with A = |0) (1| + |1) (0| + |2) (2|. 
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Then, the transformed state is p' (0) = Up (0) U'^ = 
|P+ + + where the projectors are 

P+ = |^l/+)(vl/+|,|vl/+) = -ij(|00) + |ll) + |22)),p+ - 

i(|01)(01| + |12)(12| + |20)(20|),p_ 
i(|10) (101 + 121) (211 + 102) (02|). 

Evidently, the states p' (0) and p (0) are equally use- 
ful quantum resources without decoherence. In Ref.fs^, 
Horodecki demonstrated that p' (0) is a free entangled 
state for 4 < a < 5. It will evolve into p' {t) according 
to Eq.(l), whose partial transpose always has a negative 
eigenvalue in finite time t. So p' (t) never experiences 
ESD when subjecting only to the local dephasing noise 
in Eq.(l). Furthermore, the state p' (t) does not undergo 
DSD in finite time. To see this, we perform the local op- 
eration (|1) (1| + |2) (2|)(g)(|l) (1| + |2) (2|) onp'it). The 
resulting two-qubit state is easily proved to be entangled 
and hence distillable in finite time t. So the state p' (t) 
can always be distilled under the local dephasing noise. 

By far we have considered DSD in the presence of 
local dephase noise. Another useful physical quantity 
in quantum-information problems is the fidelity, which 
measures to what extent the evolved state is close to 
the initial one. For concreteness, we study the states 
p{t) and p'{t) by means of the Bures fidelity [3§|. The 
Bures fidelity of states p and a is defined as F (p, a) = 

[^^ (V VP'^y/p)] • state p {t), the fidelity is given 

by 




while for p' (t) it reads 



6e-" ( 1 + 2e^t ^ + ge 
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(6) 



Eq.(5) and Eq.(6) indicate that the fidelities do not de- 
pend on the coefficient a. In Fig. 3, we plot F^ (t) and 
FP (t) as the function of t with F = 1. Obviously, 
the fidelity F^ (t) is always larger than F^ (t), which 
means that p (t) degrades faster than p' {t) . In the infinite 
time limit, the values of F^ (t) and F^ (t) will approach 

to [fi (15 + 2V6)]^ and (15 + 2^/3)] ^ respectively. 
Thus, we have shown that p' {t) have the advantages over 
p{t) in two aspects: Firstly, it does not experience DSD 
in finite time. Secondly, its fidelity decays slower than 
that of p {t) for all time. 



III. CONSTRUCTION OF DSD-FREE STATES 

In the last section we have presented the examples of 
qutrit-qutrit states with and without DSD. We call the 
latter as DSD-free states, which always have negative 
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FIG. 3: The fidelity versus time t for different initial states. 
The solid and dashed lines correspond to F'' (t) and F'' (t) 
with r = 1, respectively. 



partial transposition (NPT). NPT is a necessary condi- 
tion for distillability under LOCC, so the DSD-free states 
are available resource for quantum-information tasks as 
they have robust entanglement against local dephasing 
noise described in Eq.(l). It is thus worth finding out di- 
verse DSD-free states in terms of different kinds of local 
dephasing noises. In this sense, DSD-free states deserve 
storage in quantum-information "warehouse" and they 
are more useful than DSD states. On the other hand, 
one can also directly propose schemes of quantum infor- 
mation processing based on the robustness of DSD-free 
states. 

We now turn to develop several approaches to con- 
struct DSD-free states for an arbitrary qutrit-qutrit state 
a. As a direct method, we can locally project a few copies 
of the state a{t) onto some 2 (g) N{N > 2) states. Such 
states are distillable (DSD-free) under LOCC, if and only 
if they have NPT in finite time [33|. Generally, it's dif- 
ficult to find out the projectors and there have been a 
few results on this problem (33 . fiol . l4l| . Here we only 
investigate one copy of a{t) under the dephasing noise 
described in Eq.(l). For example: 

Lemma 1. Consider an NPT entangled state 
(t(0) and it experiences the local dephasing noise 
in Eq.(l). Then a{t) is always distillable in fi- 
nite time, when XlLi (0 ^1 W <^ (0) (0 ^2 {t) or 
ELi dI (t) eI (t) a (O) ^2 (t) A (i) is entangled. 

Proof. One can get the above given state by project- 
ing a{t) with the projector D2 (t) on system A or E2 {t) 
on system B respectively (both projectors actually have 
nothing to do with time t when separately used as oper- 
ators) . Since the given state is a distillable qubit-qutrit 
entangled state, the assertion follows from [33|. It's in- 
teresting to note that the given state is just the sum of 
the third and fourth terms in Eq.(l). ■ 

Lemma 1 actually requires to check the PPT of some 
6x6 matrices depending on time t, which could be done 
by using Mathematics. This is a systematic way to build 
the subspace of robust entangled states that never expe- 
rience DSD in Eq. (1). To show the connection between 
the last section and lemma 1, we propose a weak version 
of lemma 1 as follows. 

Lemma 2. Suppose {t) e\ (t) a (0) E2 {t) D2 {t) are 
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entangled, then the states fT(0) are DSD-free states in 
local dephasing noise Eq.(l). 

Proof. The state XlLi (t) sj (t) a (0) E, (t) D2 (t) is 
entangled when (t) {t) a (0) E2 (t) D2 (t) is entan- 
gled. The assertion then follows from lemma 1. ■ 

One may see that the state p'(0) in the previ- 
ous section belongs to the case of lemma 2. The 
advantage of lemma 2 is that the density operators 
{t) eI {t) a (0) E2 (t) D2 {t) do not contain time t and 
are convenient for calculation. This result is reasonable 
for the local dephasing noise described in Eq.(l) which 
occurs only between the ground state i = and the ith 
(i = 1,2) excited state. 

The model in Eq.(l) is neither the simplest case of local 
dephasing nor the most general case. In the most gen- 
eral case dephasing occurs between all local basis states 
within each subsystem. In what follows we put forward 
some results under such general local dephasing noise 

Lemma 3. Qutrit-qutrit entangled MC states are 
DSD-free states in general local dephasing noise. 

Proof. Let us consider the qutrit-qutrit maximally cor- 
related (MC) states aMc{Q), which read (i^ l 

2 

(0)= ^ a,,|«)(jj|. (7) 

We make it go through the channel of general local de- 
phasing noise. As we know, such noise cannot change the 
diagonal elements and the non-diagonal elements of the 
evolved state have non-zero dephasing values for any fi- 
nite time. So the evolved state (7Mc{t) is still a MC state. 
On the other hand, the MC state is separable if and only 
if all non-diagonal entries are zero. Then (TMc{t) is al- 
ways entangled for all finite t whenever the initial state 
fj\fc(0) is entangled. Furthermore, one may prove that 
(^Mcii) can always be distillable by using the similar ar- 
gument to show the distillability of the state p' (t) in the 
previous section. ■ 

This example enlightens us on building DSD-free en- 
tangled states in a stronger way, i.e.. 

Lemma 4. We perform the local operation 
(K) (*l + \j) 01)^(1"^) ("^1 + l*^) ('^D on any initial qutrit- 
qutrit state, where z, j, to, and n are chosen from the basis 
0, 1 and 2, respectively. If the projected state is an en- 
tangled MC state, then the initial state is a DSD-free 
state in general local dephasing noise. ■ 

One can prove lemma 4 by showing that there are al- 
ways some nonzero non-diagonal entries of the projected 
state, and we can take it as the new initial state in the 
local dephasing noise. Evidently, the newly built state in 
lemma 4 contains MC states in lemma 3. 

These results can be extended to the scenario of higher 
dimension. For simplicity, we consider the general MC 
state (j'Mci^) = Sfjio Suppose it experi- 

ences generaHzed local dephasing noise. Analogous to the 
arguments for lemma 3, one can see that the state cr'Mci^) 



also has no DSD in finite time. Thus we have provided a 
high dimensional subspace {|00), |11), , |d — 1)} 
in which all entangled states are distillable under local 
dephasing noise; that is. 

Lemma 5. Entangled MC states are DSD-free states 
in general local dephasing noise. ■ 

One can also build other higher dimensional DSD-free 
states, which could be locally projected onto entangled 
MC states by following the techniques for lemma 4. As a 
short summary, lemma 1 to 5 give a few primary results 
on DSD-free states in finite time, in terms of the spe- 
cial dephasing noise Eq.(l) and general dephasing noise 
respectively. 

To gain a better understanding of the entanglement 
dynamics, we further study the properties of the evolved 
state in the infinite-time limit, which is summarized as 
the following lemma. 

Lemma 6. For any qutrit-qutrit state under the de- 
phasing noise in Eq.(l), the final evolved state could be 
separable or eternally distillable in the infinite-time limit, 
but no PPT entangled sate can exist as time goes to in- 
finity. 

Proof. Let us briefly justify the above statement. 
Consider an arbitrary qutrit-qutrit state a under the 
action of an infinite time of local dephasing, then the 
final state a (t) can always be written as: a|00) (00| + 
|0)(0| ® (6|l)(l|+c|l)(2|+c*|2)(l|+d|2)(2|) + 
(e|l)(l|+/|l)(2|+.r|2)(l|+g|2)(2|) ® 
|0)(0| + (|1)(1| + |2)(2|) ® 

(|1) (1| + |2) (2|)p(0) (|1) (1| + |2) (2|)® (|1) (1| + |2) (2|). 
There are four terms in all, and each of them corresponds 
to one of the terms in Eq.(l), respectively. Note that 
the first three terms could be removed by the local 
operator (|1) (1| + |2) (2|) ® (|1) (1| + |2) (2|) and it does 
not change the fourth term. So the evolved state a {t) 
is entangled if and only if the fourth term is entangled. 
It immediately implies that the state a {t) could be 
separable or eternally NPT and distillable in the infinite 
time Hmit. As we have seen in the previous sections, 
both cases also exist in finite time's evolution. Since the 
entanglement property of the evolved state is determined 
by the two-qubit state, then we conclude that any PPT 
entangled state is unstable as time goes to infinity under 
the local dephasing noise in Eq.(l). In contrast, it's 
still an open question whether there always exists PPT 
entangled state in finite time's evolution. ■ 



IV. DISCUSSIONS AND CONCLUSIONS 

Before concluding we would like to discuss some open 
questions as follow: 

Firstly, the initial state in our scheme is a free entan- 
gled state. We have shown that it may evolve into a PPT 
entangled state and then a separable state for finite time; 
it may also be eternally free entangled for all finite time in 
the local dephasing noise described by Eq.(l). We have 
also shown that any qutrit-qutrit PPT entangled state 



6 



will lose its entanglement as time goes to infinity in our 
decoherence model. Then an open question is whether 
PPT entanglement can be preserved in any finite time 
under (general) local dephasing noise? If the answer is 
no (supported by the calculation in Sec. II), then we have 
another interesting difference between the free and bound 
entangled state: that is, the former can be robust against 
the decoherence in any time while the latter cannot. 

Secondly, we have proposed a few systematic methods 
to build DSD-free subspaces for 3 (g) 3 state, in which 
all states keep their distillability asymptotically in local 
dephasing noise. The primary extension to higher di- 
mensional bipartite system has also been proposed. It's 
then an open question to find out more entangled states 
which are robust against different local noises (beyond 
our model) in any dimensional space. Such states are 
applicable to quantum information schemes and thus de- 
serve a deeper study. 

In summary, we have introduced the concept of DSD 
through an explicit qutrit-qutrit state. We found that un- 
der the bi-local dephasing noise, the free entangled state 
may evolve into a PPT entangled state (DSD) and then 
become separable (ESD) permanently. Moreover, a sim- 
ple local unitary operation on the initial state can avoid 



the sudden death of distillability. We also have compared 
the action of DSD and DSD-free states in terms of the 
Bures fidelity. Next, we have proposed the systematic 
methods of building the DSD-free space against decoher- 
ence, so that the evolved states are distillable in finite 
time. Finally, we proved that there is no PPT entangled 
state in infinite-time limit in our model. Our results im- 
ply that further study on free entanglement with time 
evolution in practical local noises is required. 
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